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Abstract 

We present an explicit computation of the CP asymmetry in scattering processes 
involving the heavy right-handed neutrinos of the type I seesaw framework and the 
Standard Model gauge bosons. Compared to CP violation in two-body decays and 
in scatterings with top quarks there are new contributions at one loop in the form 
of new type of vertex corrections as well as of box diagrams. We show that their 
presence implies that, unlike the CP asymmetry in scatterings with top quarks, 
the CP asymmetry in scatterings with gauge bosons is different from the two-body 
decay asymmetry even for hierarchical right-handed neutrinos. This also holds for 
the L-conserving CP asymmetry in scatterings with U(l)y gauge bosons. 
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The discovery of neutrino oscillations makes leptogenesis a very attractive 
solution to the baryon asymmetry problem [HE]- In the standard framework 
the tiny neutrino masses are generated via the (type I) seesaw mechanism 
[3H6] and since the new singlet neutral leptons have heavy lepton number (L) 
violating Majorana masses they can produce dynamically a lepton asymmetry 
through out of equilibrium processes. Eventually, this lepton asymmetry is 
partially converted into a baryon asymmetry due to fast sphaleron processes. 

In this framework the CP asymmetry can be generated in the iV, decays [7J, 
i.e. Ni — » IjH [8J as well as in scattering processes of the iVj with particles 
in the plasma. In perturbation theory, to lowest order, the CP asymmetry in 
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any of these processes arises from the interference between the tree level and 
one-loop amplitudes. The CP violating asymmetry in scatterings with top 
quarks has been considered in Refs. [5HT2"]. In Refs. [9HTT] it was argued that 
an approximate equality between the scattering and decay asymmetries should 
hold for hierarchical right-handed neutrinos (RHNs). Ref. [12] presented an 
explicit computation of the CP asymmetry in top quarks scatterings (TQS) 
and directly showed that this was correct. It is important to note that there is 
a one to one correspondence between one-loop diagrams in decays and TQS. 
Then, as discussed in Ref. [TT], the "factorization" of a common CP asymmetry 
can be easily understood in terms of an effective field theory in which all but 
the lightest of the RHNs have been integrated out and their effect appears in a 
dim-5 operator (H£)(H£). In this approximation only "bubble-like" diagrams 
contribute to the one-loop amplitudes of both decays and TQS. 

As for the CP asymmetry in scatterings with gauge bosons, its effect was esti- 
mated in Ref. [12] under the assumption that it also factorizes in terms of the 
decay CP asymmetry. However, the same argument that leads to the under- 
standing of the equality between the CP asymmetries in decay and TQS does 
not hold for gauge boson scatterings (GBS) due to the presence of additional 
contributions to the amplitude at one loop (which in the effective field the- 
ory approximation contain three rather than two particles in the loop). With 
this motivation in this note we carry out an explicit calculation of the CP 
asymmetry in GBS . 

For a given heavy neutrino species (called here iVj) there are three different 
types of GBS, namely 

N t A tjH , Nilj HA , IjA N { H , (1) 

and each one has an associated CP conjugate process. Here A = W a or By for 
SU(2)l and U(l)y bosons, respectively. The Lagrangian for the Yukawa and 
gauge interactions relevant for the computation of these processes reads 



Y+A : 



(9 M if) t fl£# + h.c. 



+E 



(2) 



where i,j are generation indices of RHNs and lepton doublets £j = (vj,£j ) 



1 In principle one must consider all processes at a given order in the coupling 
constants. In particular the three body decays iVj —> ijHA should be included 
together with the gauge boson scatterings. However the three body decay rate is 
phase space suppressed, thus the processes iVj — > IjHA and their CP asymmetries 
can be safely neglected (see e.g. the appendix of [12J). 
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Fig. 1. Tree diagrams for the three scattering processes in Eq. ([!]). 

respectively, Pl,r = \ (1 T 7s), an d H = ir 2 H* with H T = (h + , h°) being the 
Higgs doublet. 

We denote by 7^ = j(A — > B) the thermally averaged density rate for a 
state A to go into the state B (summed over initial and final spin and gauge 
degrees of freedom), and by A 7^ = 7^ — 7^ the CP difference between the 
rates for particle and antiparticle processes. Neglecting the thermal motion 
with respect to the plasma, the relevant thermal average rates for 12-^34 
processes can be written as 




where sll(s) is the reduced cross section, K±(x) is the modified Bessel function 
of second kind of order 1, and T is the temperature. We have also introduced 
the Mandelstam variables in terms of the momenta pi of the four particles in 
the process, s — (p x + P2) 2 , t = (p± — P3) 2 , and u = (pi — p^) 2 , with p\ = mf. 
Thus for each process s_ = Max [(mi + m 2 ) 2 , (m 3 + m^) 2 ] and 



-rrin 



-m: 



As 



(s+mf — m|) 2 
4s 



■m 2 ±1 



(s+m 2 - 
47 



■m 2 ) 2 



•m; 



For the sake of simplicity in our evaluation of the CP asymmetries we neglect 
the thermal mass of the gauge bosons, and we include the thermal masses of 
leptons rri£ and the Higgs boson m H only in the propagators (but not in the 
loops), in order to regularize the infrared divergences that appear when these 
states are exchanged in a t- or u-channel. 
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In this approximation the tree level amplitudes for the three processes in (fTJ) 
are related by crossing symmetry and are given by 

\M(N t A -> i 3 H)\ 2 = \M{N t A -> I j H)\ 2 = 2g A \Y N J 2 I (s,t,u), (3) 
|M(iV^ HA)\ 2 = \M{Nd 3 -> HA)\ 2 . = -2Q A \Y Ni .\ 2 I (u,s,t), (4) 
\M(ljA iVW)| 2 = |M(^ -> AW|) 2 = 2£ A ir^f/oK^s), (5) 

where Q A = 3#f/2 or # 2 /2 for SU(2) L and U(l) y GBS, respectively 
Io(s, t, u) = 2DtnD u £sM 2 + D 2 e u [M 2 — s), and Dxp comes from the propa- 
gator of particle P in the X-channel, Dxp = {X — mp) 4 . 

The CP asymmetry in GBS arises from the interference between the tree- 
level and one-loop amplitudes in Figs. [1] and [2], respectively. In Fig. [2] we 
show as a "waving" line labeled C\, C[ or C 2 the possible branch cuts in 
which the particles in the propagators can be on-shell and therefore give a 
contribution to the imaginary part of the corresponding amplitude. As seen 
in Fig. [2] besides the Nj self-energy (wave) corrections (diagrams (a) and (b)), 
and the vertex corrections (diagrams (c) and (d)), which are also present at 
one loop in two-body iVj decays and in TQS, there are new type of vertex 
diagrams (labeled (e) and (f)), as well as contributions from boxes (diagrams 

(g) and (h)). Furthermore for process £jA — > NiH additional imaginary parts 
appear from the last two diagrams 3.(i) and 3.(j). 

For each process the sum of amplitudes of wave diagrams (a) and (b) is gauge 
invariant. The sum of amplitudes for vertex diagrams (e) and (f) is also gauge 
invariant, while the amplitudes of vertices (c), (d) and boxes (g), (h) are not 
separately gauge invariant but the sum of the four amplitudes is. Finally for 
process JjA — > NiH the sum of amplitudes 3.(i) and 3.(j) is gauge invariant. 

Thus generically the CP asymmetry from any of the three processes in (CQ) can 
be written as 

120") 1 1201) , 12 («2) 12(6) , \JjA /„n. 

e 34 +e 34 +e 34 + e 3 4 + a ^n z h, (6) 
j 

where e j^ w ' vl,v2 ' b ^ j s generated by the interference between the tree- level am- 
plitudes and the self-energy amplitudes (a) and (b), the vertex amplitudes (c) 
and (d), the vertex amplitudes (e) and (f), and the box amplitudes (g) and 

(h) , respectively. We denote by Ae%f H the extra asymmetry for the process 
JjA — > NiH from the diagrams 3.(i) and 3.(j). Notice that each of the dia- 
grams (a), (b), (e) and (f) represents two amplitudes depending on whether 
the internal lepton line is a lepton or an anti-lepton, which leads respectively 
to total L-conserving or L-violating contributions. 



12 _ 734 7 3 4__ _ 

34 "E(73 i I + 7i)" 



4 



It is straightforward to show that for any of the processes in ([I]) we have 



N t A( w )_ Nil 



(to) 



iH 



~~ M 2 ' ' fc 



,JVj ( u ')_ 



is 



-Im 



Im 



(7) 



ik 



In the above, e^ H ^ is the contribution to the CP asymmetry in Ni — > £jH 
from the Ni self-energy diagrams [8] and the first (second) term arises from 
the L-violating (L-conserving) diagrams. 

The contributions to the CP asymmetries from the vertex (labeled with su- 
perscripts (v 1), (t>2)), box (superscript(6]), and the extra piece from diagrams 
3.(i) and 3.(j) (superscript (A)) readLULil 



\M(12 ->■ 34) | 2 - \M(12 ->■ 34) | 5 



(vl),(i;2),(6),(A) 



g A a k M k Mi 



t34 ' (5, t, u) + /' 3 T (s, t,u)+9(s- Mi) (s, t, u) 



F A C jk M k Mi 



2n 



(v2) 



,/12(u2), 



X 



,(«2) 




( Klf" (s, t, u) + K'H ' (s, t,u))+ g A 6(s - M 2 k )Jll^\s, t, u) 



where / 



12(«1).W, 
34 1 



s, t, u) = Igf (s, t, u) + igf (s, t, u) comprises the contribu- 



tion from the C\ cuts of the vertex diagrams (c) and (d) (given by I\i (s, t, «)) 

and box diagrams (g) and (h) (given by I^ b \s,t,u)). I'^ b \s,t,u) con- 
tains the contribution from the C[ cuts present only in the box diagrams 
of NiA — > £jH. J^f ^ (s, u) has the contribution from the G% cuts, which 
are kinematically allowed when s > M k and they are possible only in the pro- 
cesses Nilj — > HA and £jA — > NiH. All these contributions are L-violating 
and they are present for both scatterings with SU(2) i and U(l)y gauge bosons. 



2 When evaluating £34 ^ and £34 one must specify a gauge since, as mentioned 
above, the one-loop vertex (vl) and box amplitudes are not separately gauge- 
invariant. In what follows we will work in the unitary gauge. 

3 Since the resulting expressions are somewhat cumbersome, to double check our 
results we have computed the asymmetries both by explicit evaluation of the cor- 
responding loop integrals as well as by using the Cutkosky rules that give directly 
the absorptive part of the Feynman diagrams. 
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H, H" 




Fig. 2. One-loop diagrams for the three scattering processes in ([I]). We mark with a 
waving line the possible branch cuts which contribute to the imaginary part of the 
corresponding amplitude. 

J34 (s,t,u) has the - L-conservmg - contribution from the C 2 cuts in the 
diagrams 3.(i) and 3.(j) present only for IjA — >■ NiH. 

Kll^ v2 \s,t,u) (K^ (s,t,u)) contains the contribution from the C\ cuts of 
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the L- violating (L-conserving) vertex diagrams in graphs (e) and (f). Similarly, 



K'l^ v2 \s,t,u) (K'sf (s,t,u)) involves the contribution from the C[ cuts of 
the L-violating (L-conserving) diagrams in the graphs (e) and (f) of process 
NiA — > £jH. Notice that the diagrams (e) and (f) correspond to amplitudes 
where the gauge boson is coupled to the one-loop self energy Ni — N k . Since 
the RHNs are SU(2) L singlets, the sum over the SXJ(2) L degrees of freedom 
cancel in each of these graphs. This, however, is not the case for the U(l)y 
gauge boson. Correspondingly the coupling factors Ta — for A = W a while 
J 7 a = Ga = 9l/2 for A = B Y - 

The C\ cuts are kinematically allowed for all processes and their corresponding 
contributions to the amplitudes are related by crossing symmetry. For the 

amplitudes of diagrams (c) and (d) this reads lf i ^ Vl (s,t,u) = I^(s,t,u), 
Iah {V1 (s,t,u) = -lW(u,s,t), and 1%^' (s,t,u) = I^ vl \u,t,s), with 



(«2) 



I {vl) (s,t,u) 



(9) 



D 2 M 2 



su 2 



1-u 



,~A2 



a k 



+D tH D ul M 2 



1 - 



1-u 
a k + 1 — i 



tit 



L 2 (t) 



1-u) 2 
\ s 



a k + 1 - u 
1 : = L 2 (U 



1 



1-U 

a k + 1 — u 



U(u. 



l-t\ 1-t " y ') l-u\ 1-u 

The C2 contributions to the vertex asymmetry from diagrams (c) and (d) are 



[s,t,u) = D M „DttM. 



u 



1-s 



(vl) 



J^; H (s,t,u) = D tH D se M< 



u 



1 - s 



a k a k + 1 - s 

1 - — : —L A (s) 

s 1 — s 

„ a k + 1 - s 

a k — s -\ — L 4 (s) 

1 — s 



+D 2 ,M 2 



'1 - § 



S\2 



,7 ~m ~\ , t(a k + 1 - s) - a k su 

(t - u)(a k - s) H L 4 (s) 

1 — s 



(10) 



and J, 



NiA( vl ), 



s,t,u) = 0. In Eqs.(l9l)- (fT0l and the following we define s 



Mf ' 



t = M?> S = alld 



a k + x a k + 1 - x x(a k + l—x) 

Li(x)=m , L2\x) = In , L4(a;)=ln — 



x 



a k 



a k 



T ( xy + a k {l-x) xy + a k (l-x) 
L${x,y) = In — , L 5 (x,y) = In . (Ill 



a k (l - x) 



For the C% contributions to the asymmetry from diagrams (e) and (f) the 
crossing; symmetry reads Kp£ (v '(s,t,u) = K (v2 \s,t,u), K 



AH 



(v2) 



,S,t,U) 
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-K (v ^(u,s,t), K£ H (s,t,u) = K^(u,t,s), and equivalent^ for the K 
functions, with 



K^ v2 \s,t,u) = - 



2s [D tH t [(t - 2u) + u(l - s)] - D U (U [si + u(l - 5)1 } 



(s - a k )(t + u 

_ 2su {D tH t [3 - 2g + g] + -P-ui [(2t~ - j- g + u 2 ] } 
ii (s, t, u) — ~ _ . (12) 



(s - a fc )(£ + w) 3 



The contribution from the C[ cuts in diagrams (e) and (f) (which are only 
present for NiA — > £jH) is 



K'^ffs, t, u) = -K^ (s, t, u) , Kf£ ' (s, t, u) = -K^ (s, t, u) , (13) 



■<v2) 



and hence the contributions from vertices (e) and (f) to N^A — > £jH cancel 
exactly. 

For the box diagrams: I^^ b \s,t,u) = I^ b \s,t,u) = I bl (s,t,u) + I b2 (s,t,u), 
lAH iib \s,t,u) = -I {b \u,s,t), and (s,t,u) = I<V(u,t,s), with 



hi (s, t, u) = 



+D ui 



(1-5) 2 u(l-f) 

st a k + s 



Qfc + g a fc (afc + 1 - t) 
■u si + a fc (l — s) 



L 3 (s,t) 



-1 - 



(1-5) 2 l-t 



J 62 (s, u) 
D tH 



St OLu + 1 — U 

+ " 



(1 - sf 

{ St 



t 

u 



st + a k (l - s) 



sL 2 (u) (1 - s){a k + s) + a k t 

— H rz— 7; — ^3(,s, WJ 



(14) 



1 — u su + a k (l — s 

(a k + 1 — u) st — u(l — u) 



id 



(1 - S) 2 t(l - s) 



t 2 (l-u) 



1 - 3s - 2a k ^ (l-s)(a k + s) | a k - s a k (a k + s) 



t 



t 2 



s s 2 u + a k s(l — s 



L 2 (u) (15) 
L 3 (s,u) \ . 



The corresponding contributions from C[ cuts read I'^'h ( s > ^ u ) = 

!'%% (s,t,u) = and/'^ (6) (s,t,w) = I' {b) (s,t,u) = I' bl (s,t,u) + I' b2 {s,t,u), 
with 



(b) 
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D 



tH 



su a k + s 

H = — LAs) - 



I'bi (s, t, u 
D tH ' 



(i - sy 

si 



u 



dk + s a k (ak + 1 - t) 
u si + a k (l — 5) 



L 3 (s,t) 



(1-5 



Z\2 



ak + s {a k + s)(a k + l-t) ~ 
— ~ — L n s ) ~r~i 7i ^ ^sySjt) 



st a k + 1 — u 
+ — — ~ 



Z\2 



(1-5 



[ s 1 — s + st 
+D ui \ -~ - — — 4 



st + a k {l - s) 



i , - , (1 - s) (a k + 5) + a k t 

L li S ) rr.- , - M =vi L A S -, U ) 



(16) 



[su + a k (l — §)] 



(17) 



t (1-5 



Z\2 



(l-u) 2 (a k + s) 1 — u a k + l-u 



si 2 



+ 



t 



f: 2 



+ 



1 - 3s - 2ctfc (1 - s)(a k + s) a k - s a k (a k + s) 



t 



t 2 



s 2 u + a k s(l — 5) 



Li(S) 
L 3 (s,u) 



From the Eqs. (ll4p - (fIT]) . we see that for the scattering process NiA — > £jH 
the contributions from C\ and C[ partially cancel each other 0- 



The C2 contributions to the box asymmetry are J^H ib \s,t,u) = and 



J AH (S,t,u) = 

u(a k + 1-5) 



D 



sH 



t(s-l) 



LAS) + 



a k + u Ok{a k + l-s) 



t su + a k {l — u) 



L 5 (u,s) 



(18) 



+D 



(l-t)(a k -s) (a k + u) (a k + 1 - s)(a k + u) 

-L 4 {S) H zz—, — —7^ — L 5 {U,S) 



s 2 



5-1 



su + a k (l — u) 



j n% { s i t, u) = D si 



s - a k (a k -s)(a k + u) 

r^— — y zt-L 6 {u, s) 

t su + a k {l — u) 



a k + 1 - 5 



t 2 



tu — 5(1 — s) 
1^5 



L 4 (5) + (5 - t)L 5 (u, s) 



(19) 



+D 



tli 



u(a k + l-s) ._, (flfe + 1-5) (a k + u)(l - u) + a k t 
■LAs) 



i(l — s) ' t [su + a k (l — u)] 

Finally the C2 contributions from diagrams 3.(i) and 3.(j) read 



LAS) 



Jm A u (s,t,u) = -D 



M 2 (~s-a k ) 



NiH b i U J — ^st 2g 

x [D tH u + D se s(l - u)} . 



s - a k 



+ D se M 2 (~s + a k ) 



(20) 



4 By definition I 3 2 , I'W , and j|| must be real. It can be easily verified that 
in their expressions the imaginary contributions from negative arguments in the 
logarithms in Eq. (fTT|) always cancel out. 
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Fig. 3. Ratio between the L-violating asymmetry in scatterings e s = €34 for the 
different processes as labeled in the figure and the corresponding quantity for de- 
cays e D = £^ l H - To illustrate we have used mi = 0.18T and run = 0.34T. The 
left panels show the ratio as a function of T/M\ for three different hierarchies 
M2/M1 = 2,10,100. The right panels show the ratio as a function of M2/M1 for 
two values of T/M\ = 1, 5. 

The ratios between the L-violating asymmetries in GBS and the L-violating 
asymmetry in decays are shown in Fig. [3] as a function of T jM\ (left panels) 
and of M 2 /Mi (right panels). For simplicity we have considered the virtual 
effects of only one neutrino species = N 2 different from iVj = Ni. It is 
apparent that for Niij — > HA and IjA — > N\H this ratio deviates from 
unity even for hierarchical neutrinos lj. In the relevant temperature range 

5 Unlike the CP asymmetry in decays, the CP asymmetry in scatterings depends 
on the temperature, thus when comparing both types of asymmetries at a given T, 
the RHNs are considered to be hierarchical enough if 3> Mi,T. 
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for leptogenesis 0.1 < TjM\ < 10 the deviation can be of several tens of 
percent, e.g. 50% for £jB Y ->■ N { H at T = M x . Conversely for N X A ->• ijH 
the ratio tends to one for hierarchical neutrinos, hence the asymmetry can be 
well approximated by the decay one if M 2 /Mi 3> 10, but corrections appear 
at high temperatures for milder hierarchies. 

These results can be easily understood analytically by expanding to lowest 
order in Mj/M^ the expressions in Eqs. (Q, (fl2l) . and (|T4|) - (fl7|) which gives 



I (vl) (s,t,u) 



I {b) (s,t, 



u 



1 1 (s,t,u) 
2 Ml 

Mf s [D u£ u 



3t + u(l-s] 



D tH t 



t - 2u )+u(l 



Ml 



(t + u) 3 



r ib \s,t,u) = -i {b \s,t, u ), 

K^ 2 \s } t,u) = 2I^(s,t,u) 



-K'W(s,L 



K^ 2 \s,t, 



u 



Mf 2Su [D tH t [3 - 2(f + u)] + [(2t - u) - t 2 + u 2 ] } 



Ml 



(t + u 



i~A3 



(21) 

,(22) 

(23) 
(24) 

,(25) 



Including these results in Eq.flE]) and using Eqs. (J3j) — (EJ) one finds that 



NiA( vl ) 



>1) 



z N t H 



(vl) 



l^[Y N Yl 



M k 



-ijH ■ 



(26) 



In the last equality we have used that the vertex CP asymmetry from the 
decay iVj ->■ ijH 



c Ni 0) 
-ijH 



k^i 



87T (y n y£). 

\ / 3 



1 - (l + a fc )ln 



1 + Qfc 



>r- r Mi 
^ jk M~ 

k=ii ™ 
167T 



(27) 



where the last approximation holds for ^> 1. So, as in the case of 

TQS, both wave and vertex contributions to the CP asymmetry in scatterings 
with SU(2) L gauge bosons are equal to the corresponding contributions to the 
decay CP asymmetry in the hierarchical limit (M^/Mj 3> 1). 

For the process NiA — > ijH Eqs. (!22|) - (|24|) imply that in the hierarchical limit 
the CP symmetries from vertices (e) and (f) with U(l)y gauge bosons, and 
from boxes with either SU(2)x or U(l)y gauge bosons vanish, 



r NiW a ( b ) _ MiB Y ( b ) _ NiB Y ( v2 ) 



0. 



(28) 
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However for the other two processes we find that the L- violating contributions 
to the CP asymmetries satisfy 




^0, 




e jBY ( b ) _ 1 ijBrW 
e N,H - 2 e N t H 



^0. 



(29) 



We note that the vanishing of ef^ and e; 



at leading order in Mi/Mk 



is due to the presence of the additional C[ cuts, which, as mentioned above, 
results in a partial or total cancellation of the contributions from the C\ cuts 
so that only terms of higher order in Mi/Mk remain. 

We have mentioned before that for each process the sum of amplitudes of wave 
diagrams is independently gauge invariant, as it is the sum of the amplitudes 
for the vertex (v2) diagrams, but the amplitudes of vertex (vl) and box di- 
agrams are not separately gauge invariant. However at the leading order in 
Mi/Mk, the box and vertex (vl) diagrams are independently gauge invariant, 
thus the results in Eq. (l26|) and Eq. (l29|) are also gauge independent. 

Summarizing, we find that even for very hierarchical heavy RHNs 



where e t - H = e^ H + e^ H . Hence, we conclude that the CP asymmetry in 
the scattering NiA — > £jH factorizes in the hierarchical limit. However, this is 
not the case for the scatterings Nilj — > HA and £jA — > NH. 

Also, as mentioned above, in the processes Nilj — > HBy and IjBy —> NiH 
there are new L-conserving vertex contributions to the CP asymmetry. They 
imply that in the hierarchical limit the L-conserving CP asymmetry in these 
scatterings is not equal to the L-conserving CP asymmetry in decays. Note 
that although the L-conserving asymmetries are always suppressed by a higher 
power of Mi/Mk compared to the L- violating ones, they are the only source 
of CP violation in models with conservation of L [T3HT5] . 

Finally, a few comments about the impact of our results for the determination 
of the baryon asymmetry are in order. Although the CP asymmetry in GBS has 
more interesting features compared to the one in TQS, we have found that still 
the CP asymmetry per scattering is (0.5 — 2)x the CP asymmetry per decay. 
Therefore the generation of CP asymmetry through GBS is relevant at high 
temperatures (T > Mi) when the scattering rates are larger than the decay 
rates. Hence it can have effects in models for leptogenesis that are sensible to 
the high temperature regime. This includes standard type I leptogenesis in the 
weak washout regime and low energy models in which the sphaleron processes 



t e j H — t ljH T t AH T t 



(30) 
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freeze out before the RHNs have decayed [TH[T6] . 
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